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ABSTRACT 

Context. We report on XMM-Newton observations of SN 1987 A in the Large Magellanic Cloud. 

Aims. The large collecting area telescopes together with the European Photon Imaging Cameras (EPIC) provide X-ray spectra with 
unprecedented statistical quality and make it possible to investigate the spectral evolution during the brightening observed since the 
discovery in X-rays. High resolution spectra from the Reflection Grating Spectrometers yield a complementary view and allow us to 
perform more detailed investigations of prominent emission lines. 

Methods. The X-ray spectra were modeled with two-temperature emission components from a hot plasma in coUisional ionization 
equilibrium and in non-equilibrium (NEI). 

Results. We find a temperature for the equilibrium component of 0.24±0.02 keV in January 2000 and April 2001 which increased 
to 0.30±0.02 keV in May 2003 and also an indication for a temperature increase in the hot NEI component from ~2 keV to ~3 keV. 
Emission line ratios inferred from the RGS spectra suggest temperatures as low as 100 eV and an increase in the ionization state 
of oxygen and neon consistent with the observed temperature increases. The fast readout of the EPIC-pn instrument yields X-ray 
fluxes free of CCD pile-up effects which we used to normalize pile-up corrections for the published Chandra fluxes. The corrected 
X-ray light curve of SN 1987 A in the 0.5-2.0 keV energy band is best represented by a linear increase up to about day 4000 after the 
explosion and an exponential rise afterwards until the last published Chandra observation on day 6716. Modeling the light curve by 
emission from the inner ring which is approximated by a circular torus a central density nn = l.lSxlO"* cm"^ is found. In this model 
the forward shock has just passed the center of the torus. 

Conclusions. SN 1987 A continues to brighten exponentially in soft X-rays. The X-ray spectra can be represented by pure thermal 
emission without significant contribution from a compact object yet. 

Key words. ISM: supernova remnants - supemovae: general - supernovae: individual (SN 1987 A) - X-rays: general - X-rays: stars 



1. Introduction 

The first supernova observed in 1987 (SN 1987 A) went off in the 
Large Magellanic Cloud (LMC) on 23 February. This event with 
observations involving the entire electromagnetic spectrum and 
the first detection of cosmic neutrinos brought about a number of 
astrophysical surprises, some of which are puzzling still today. 
In this paper we concentrate on the observations of soft (0.1 - 3 
keV) X-ray emission since the explosion, thought to arise from 
the interaction of shocks with the ambient circumstellar medium, 
and possibly a central compact object like a neutron star We 
present and summarize for the first time the full set of previous 
data, which to some extent have been re-calibrated, and we add 
the XMM-Newton observations. 

The first attempts to search for soft X-rays were made in 
August 1987 with sounding rockets caiTying imaging telescopes. 
An upper limit of 1.5x10^^ erg s ' could be established which 
excluded an ambient dense and slow wind typical for a red 
supergiant wind, but which is still compatible with a blue su- 
pergiant ( Aschenbach et al. 1987). X-rays from SN 1987 A be- 
came clearly visible in the ROSAT images not much before 
spring 1992 ( Beuermann et al. 1994; Gorenstein et al. 1994). 
The monitoring of SN 1987 A with ROSAT revealed a slow, but 
steady increase in flux over the first years ( Hasinger et al. 1996). 
The brightening continues since theri but with inc reasing pace. 
iBurrows et al liiOOtlPark et alJ2002il2()04ll20()5h . Chandra has 
taken images of SN 1987 A since late 1999, which show a ring- 
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like emission region steadily growing in diameter. Details of the 
ring structure became apparent in the meantime. These structural 
features are coiTelated with the optical hot spots, first discov- 
ered in earlier images taken with the Hubble Space Telescope 
(HST). These optically bright spots evolve quite rapidly and 
they now seem to completely engulf the so-called inner ring 
jLawrence etalJ20 00: McCrav 2005). 

The ring with its fairly high matter density existed before 
the explosion and was detected by re-emitted light of the prompt 
supernova EUV flash. A substantial heating, up to X-ray tem- 
peratures, of the ring was predicted for the period between 2002 
and 2010 when the supernova shock wave would have reached 
and entered the ring. Both the optical and X-ray images and the 
X-ray light curve demonstrate that this state has been arrived at. 
Hope is that we can learn more about the details of the interac- 
tion process, the origin of the ring as well as of the progenitor 
star and its evolution by X-ray spectroscopy. 

The medium energy resolution X-ray spectrum measured 
by the Chandra ACIS detector is composed of two emission 
components with temperatures kT of ~2.4 keV and ~0.22 
keV, which sometimes are interpreted as contributions from the 
blast wave s hock and the decelerated slow shock, respectively 
(Mi chael et al . 2002; Park et al. 2004). The strong brightness in- 
crease made it also possible to obtain high resolution X-ray spec- 
tra with the Chandra LETG. A total exposure of 289 ks was ac- 
cumulated between August 26 and September 5, 2004. InfeiTed 
temperatures from the emission line ratios range from 0.1 to 2 
keV an d expansion velociti es are found between 300 and 1700 
km s-i (IZhekovet alJ200.5h . 
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Table 1. XMM-Newton observations of SN 1987 A and data selection. 



Instrument 


Read-out 


Filter 


Sat. Re vol. 


Date 


Time (UT) 


Net Exp.' 


Counts^ 




Mode 










[s] 






FF, / J ms 


JVleQitim 


9 1 


ZUUU Jan. Ly/ Z.V 


nA-9S 




9 1 OA 

z. Ly^ 




FF, / J lliS 


iVlCUiUlll 


99 


9nnf) Tan 91 
iivvuu J an. z. i 


i J. J>0 — 1 7. JT- 




J i 


EPIC-pn 


FF, 73 ms 


Thin 


22 


2000 Jan. 21/22 


20:32- 12:02 


30990 


3060 


EPIC-pn 


FF, 73 ms 


Medium 


244 


2001 Apr. 8/9 


20:35-07:14 


24534 


3952 


EPIC-pn 


SW, 6 ms 


Medium 


626 


2003 May 10/11 


16:41-23:37 


55855 


25609 


EPIC-MOSl 


FF, 2.6 s 


Medium 






16:35-23:37 


101520 


16020 


EPIC-M0S2 


FF, 2.6 s 


Medium 






16:35-23:37 


101510 


16018 


RGSl 


Spectro 








11:42-23:38 


112160 


2509 


RGS2 


Spectro 








10:47-23:38 


111670 


3358 



'Net exposures for the observations from revolutions 21/22 were derived from a comparison of the spectrum of N 157 B with that obtained from 
the ROSAT PSPC due to missing timing information in the early data (Haberl et al. 2001). 

^Net counts as used for the spectral analysis in the following energy bands: 0.2-9.0 keV for the EPIC instruments (except for EPIC-pn during 
revolutions 21/22 when we adopted 0.3-9.0 keV due to the higher low-energy threshold used at the beginning of the performance verification 
phase); 0.3-2.0 keV for RGS. 



One other outstanding issue is the existence and the time 
of appearance of a central compact remnant like a neutron star 
So far the results have been negative. The most stringent up- 
per limits for an X-ray point source in the newly born super- 
nova remnant are based on the analysis of Chandra images in 
the 2-10 keV band with < 1.5x10^^'* erg s"' (Park et al. 2004) 
and on INTEGRAL measurements in the 20 -60 keV band (Lx 
< LlxlO^*" erg s IShtvkovskiv et alJl2005^. The H ST derived 
optical hmit is Lopt < 5x10^^ erg s~' ( iGraves et alJ2 005). 

XMM-Newton observed SN1987A on several occasions. 
Here, we report the results of the spectral analysis of the data 
with the highest statistical quality available to date, making use 
of the EPIC cameras and the RGS spectrometers. 

2. X-ray observations 

XMM-Newton iTansen et al.l 1200 ll) observations from three 
epochs were used to investigate the spectral evolution of 
SN 1987 A (for observation details see TableQ. The 30 Doradus 
region in the Large Magellanic Cloud was selected for first- 
light including SN1987A at an off-axis angle of about 7 arc 
minutes (Dennerl et al. 2001). Five e xposur es were taken with 
the EPIC-pn instrument ( Struder et alJ l200ll) in full-frame (FF) 
imaging mode which we use for spectral analysis as described 
by p^aberl et al. (2001). Exposures with the same instrumental 
setup (pointing direction, optical blocking filter) were merged 
resulting in three EPIC-pn "first-light" spectra of SN1987A 
from January 2000. An observation at a second epoch in 
April 2001 was performed utilizing all instruments on board 
of XMM-Newton ; the CCD cameras EPIC-pn and EPIC-MOS 
llTurn er et al. 2001) and the Reflection Grating Spectrometers 
(RGS, iden Herder et aL.2001.) . For spectral analysis we use the 
EPIC-pn data from this observation. 

Spectra from all instruments were accumulated from the 
~100 ks observation performed in May 2003 when the source 
was sufficiently bright such that the long exposure yielded 
spectra with sufficient statistical quality for all instruments. 
During this observation EPIC-pn was set to the faster small- 
window (SW) CCD readout mode. For optical light blocking the 
"medium" filter was selected in the EPIC cameras, except during 
the last part of the first-light observations, when a thin filter was 
used. 



We screened the data for high background intervals and ex- 
tracted source and background spectra of SN 1987 A from cir- 
cular regions with a 25" radius placed on the source and a 
nearby point-source free area. For the EPIC-pn spectra single- 
pixel (PATTERN 0) events were selected while for EPIC-MOS 
all valid event patterns (PATTERN to 12) were used. The spec- 
tra were binned to contain a minimum of 30 and 40 counts per 
bin for EPIC and RGS, respectively. Table[2Hsts net exposures 
and net counts accumulated for the spectra. The data from satel- 
lite revolutions 244 and 626 were processed with SAS release 
6.1.0 together with the latest calibration files relevant for spec- 
tral analysis of EPIC-pn data (EPN_REDIST_0010.CCF for the 
re-distribution matrices and XRT3_XAREA_0010.CCF with an 
update to the mirror eff'ective area). The spectr a from revolu- 
tions 2 1/22 were directly taken from the work of 'Denne rl et alJ 
(1200 I h and Haberl et al. (2001) together with new response files 
mentioned above. Spectral fitting was performed using XSPEC 
(iArnaud»1996) version 11.3.2i. 



3. Spectral analysis 

For studying the spectral evolution of SN1987A we concen- 
trate on the EPIC-pn data because they provide the spectra with 
the highest statistical quality, in particular for the early obser- 
vations when the source was still comparatively faint. We si- 
multaneously fit the five EPIC-pn spectra with different two- 
temperatu re models similar to those applied to Chandra ACIS 
spectra bv'Pa rk et al.l (l2004ft . The best results were found using 
a two component thermal model (c.f. Table |5J involving both 
the collisional ionization equilibrium models (XSPEC models 
MEK AL or RAYMOND, M ewe et al., 1985.; Raymond & Smitll 
|192|) and/or the non-equilibrium model (NEI, iBorkows ki et alJ 
'2001"). For the photo-electric foreground absorption we use a 
fixed Galactic column density of Nh - 6x10^" cm"^ with the 
elemental abundances given by Wilms et al. ( 2000) and an LMC 
column density (free in the fit) with the same abundances rela- 
tive to H but reduced by a factor of 2 for each element except He 
(also relative to the Wilms et al. abundances). The, in compar- 
ison with the solar value of Anders & Grevesse (J.989), lower 
Oxygen abundance of Wilms et al. better fits the EPIC spectra 
around the Oxygen K absorption edge. 
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Table 2. Fit quality of two-temperature models. 



Model' 






EPIC-pn, all 


MOS 


, rev. 626 


RGS 


, rev. 626 










dof 


9 

X'r 


X' 


dof 


x; 


x' 


dof 


9 

Xl- 


A: PHABS 


* VPHABS 


* (VNEI + VRAYMOND) 


806.5 


716 


1.126 


864.1 


357 


2.434 


799.0 


580 


1.377 


B: PHABS 


* VPHABS 


* (VNEI + VMEKAL) 


804.8 


716 


1.124 


903.6 


357 


2.545 


811.4 


580 


1.399 


C: PHABS 


* VPHABS 


* (VNEI + VNEI) 


802.2 


715 


1.122 


890.4 


357 


2.508 


870.2 


580 


1.500 


D: PHABS 


* VPHABS 


* (VMEKAL + VMEKAL) 


858.1 


717 


1.197 


878.4 


355 


2.474 


810.3 


580 


1.397 



'Model components as available in XSPEC: PHABS represents the Galactic foreground absorption, VPHABS the absorption in the LMC with 
reduced metal abundances, VRAYMOND and VMEKAL emission from hot plasma in collisional ionization equilibrium with variable abundances 
and VNEI emission from a non-equilibrium plasma (see also Sect. 3). 




Fig.l. EPIC-pn spectra of SN1987A 
from Jan. 2000, April 2001 and May 2003 
(in increasing intensity). The histograms 
show the best fit model A as defined in 



Channel energy (keV) Table and parameters listed in Tables 

andH 



We allowed some key parameters of the models to vary with 
time, i.e. they are free parameters for the individual spectra 
of different epochs. These are temperature and normalization 
(given in the following as emission measure). The ionization 
timescale in the NEI models was found to be the same within the 
error limits for the various spectra and therefore it was used just 
as one single free parameter in the combined fit. Similarly, we 
assumed that the set of elemental abundances used in the emis- 
sion models did not change with time as a whole, but we still 
fit abundances individually for each element. Also for the emis- 
sion models abundances relative to those given by Wilms et al.l 
( 1200 0) are used. For the three first-light spectra only a constant 
factor was allowed to vary to account for the uncertainties in the 
net exposure time. However the spectra agree within ±3%. 

The energy resolution of the EPIC instruments does not al- 
low to differentiate between the different two-temperature mod- 
els. Table 12 summarizes the quality of the fits for the best four 
models. The combined EPIC-pn fits yield a somewhat smaller 
reduced values for models A, B and C which each involve 
a non-equilibrium model for the high energy part of the spec- 
trum. The EPIC-pn spectrum together with the best fit model A 
is shown in Fig. [2 We cross-checked the results from the EPIC- 
pn fits with the EPIC-MOS and RGS spectra from revolution 626 
(earlier spectra have insuflicient statistics) by folding the model 
through the corresponding detector response (S AS 6.1.0 version) 
but optimizing the normalization through a constant factor The 



MOS spectra (Fig. |5} are fully compatible with the pn results 
above ~0.8 keV. Below that energy the MOS fits suffer from the 
changing event redistribution which is not included in the cal- 
ibration files of SAS 6.1.0. The RGS spectra with their higher 
energy resolution (Fig.O formally exclude model C (this model 
produces a strong line at 0.512 keV which does not exist in the 
data). 

In the following we present more detailed results for models 
A and B. The best fit parameters are summarized in Tablesl^and 
13 Model A is a combination of NEI and RAYMOND emission 
components while model B uses NEI and MEKAL. The medium 
resolution EPIC-pn spectra can not discriminate between these 
two models, but we prefer RAYMOND rather than MEKAL be- 
cause the latter produces a S XIV line at 0.512 keV which is not 
seen in the data and not present in any other model (see RGS 
spectra in Fig.O. Apart from this difference the best fit param- 
eters of model A and model B agree very well within their error 
limits. 

For both models there is clear evidence that both the tem- 
perature of the soft (MEKAL/RAYMOND) component as well 
as the temperature of the hard (NEI) component have increased 
after the second epoch (rev. 244) observation. This can also be 
directly seen in the shape of the EPIC-pn spectra below 1 keV 
where line intensity ratios differ between the spectra from rev- 
olutions 244 and 626 (Fig.[0. In particular the NelX triplet in- 
creased in strength with respect to O VIII. The increase in tem- 
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0626 EPIC-MOSl FF medium 
0626 EPIC-M0S2 FF medium 




0.5 



1 2 
Cliannel Energy (IceV) 



0626 RGE 1 
0626 RGE 3 




0.5 



Fig. 2. EPIC-MOS spectra of SN 1987 A 
from satellite revolution 626. The his- 
tograms show the prediction of model A 
with the best-fit parameter values of the 
EPlC-pn spectrum except the normaliza- 
tion (a constant scaling factor of 0.98 and 
0.99 for MOSl and M0S2, respectively). 



Channel Energy (keV) 



Fig. 3. RGS spectra of SN1987A from 
satellite revolution 626. The model is as for 
Fig|2with scaling factors 1.01 and 0.98 for 
RGSl and RGS2, respectively. 



perature may be related to the progressing equilibration of elec- 
tron and ion temperatures. 

The EPlC-pn derived best fit model A reproduces most of 
the emission lines resolved in the high resolution RGS spec- 
tra, except the O Vll He-like triplet near 574 eV (Fig.|3}, which 
suggests that some plasma with temperatures lower than 0.3 
keV contributes to the emission. Oxygen line ratios derived 
from the September 2004 Chandra LETG spectrum, taken about 
1.5 years after the last XMM-Newton observatio n, are consis- 
tent w ith lower temperatures of about 100 eV dZhekov et al.l 
l2005h . For a quantitative comparison of the line strengths of 
the most prominent emission lines seen in the LETG and 
RGS spectra we modeled the RGS spectra in a way simi- 
lar to the work of Zhekov et al. (2003) but with the exception 
that we modeled both RGS spectra simultaneously with an ab- 
sorbed bremsstrahlung component for the continuum and a set 
of emission lines with Gaussian line profiles. For complete- 
ness, we included all those emission lines which are detected 



in the RGS spectra of the bright Magellanic Cloud supernova 
remnants 1E0102.2-7219 (Rasmussen et al. 2001) and N132D 
(Behar et al. 2001). This makes sure that no observed emission 
line is missed in the model as is the case for the plasma models 
used above. The width of all lines was forced to be the same, i.e. 
one free parameter. This resulted in a best fit value for the full 
width at half maximum of 5.3±1.0 eV which corresponds to a 
Doppler velocity of 2430+460 km s ' . 

All line energies were fixed at the laboratory values, but red- 
shifted by the LMC Doppler shift of 286 km s ', which was 
found to be consistent with the LETG spectrum (Zhekov et aO 
2005). Leaving the common red-shift for all lines as free param- 
eter in the fits a Doppler velocity of 460+126 km s"' (90% confi- 
dence) is inferred from the RGS spectra. Taking into account the 
RGS wavelength scale calibration uncertainties of ~200 km s"' 
the LETG and RGS derived LMC Doppler shifts are consistent 
with each other The line fluxes derived for the fixed Doppler 
velocity of 286 km s"' are listed in Table |5] We assumed 286 
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Table 3. Results from the two-temperature model fits to the 
EPIC-pn spectra - the parameters allowed to vary with time. 
Fluxes are given in units of 10"'^ erg cm"^ s"' and emission 
measures EM in 10'''^ cm"-'. 



Parameter 


Rev. 626 
May 2003 


Rev 244 
Apr. 2001 


Rev. 21/22 
Jan. 2000 


Model A 
NEI kT (keV) 

NEI EM 

RAYMOND kT (keV) 

RAYMOND EM 

Flux (0.5-2.0 keV) 
Flux (3.0- 10.0 keV) 
Flux (0.2-10.0 keV) 


3 99+0.53 
•'•^ -0.49 

1 00+°''* 

^•""-0.15 

30+° ''2 

"■-^"-0.02 

8.10 

I. 77 

II. 4 


1 76+°-6** 

69+° '** 

24+°-"- 

5 1+'-' 

2.77 
0.43 
3.82 


2 10+°" 

-0.29 

43+° °* 
24+° °2 

3 9+0.7 
-'•^-0.9 

1.75 
0.39 
2.59 


Model B 
NEI kT (keV) 

NEI EM 

MEKAL kT (keV) 

MEKAL EM 

Flux (0.5-2.0 keV) 
Flux (3.0- 10.0 keV) 
Flux (0.2- 10.0 keV) 


-1 1 1 +0.63 
-'•^'^-0.45 

1 03+°''' 

Q3J+0.02 

8.11 

I. 74 

II. 4 


1 65+°"- 
74+°-25 

"• '^-0.15 

25+°-°^'' 

4 0+°' 

2.78 
0.40 
3.79 


2 01+017 
45+° °'' 

"•^■^-0.05 

24+° °' 

"•^^-0.02 

1.75 
0.38 
2.58 



Table 4. Results from the two-temperature model fits to the 
EPIC-pn spectra - the parameters assumed to be constant with 
time. 



Parameter 


Model 
A 


Model 
B 


Galactic Nh (10^' cmT^, fixed) 
LMCNh (10-' cm-2) 

net(I0'° scm-'') 


0.6 

2 50+°" 

^■-'"-0.33 

5 7+' "' 


0.6 

9 97+0.19 
^■^'-0.22 

5.8!li 


He (fixed) 


2.57 


2.57 


C (fixed) 
N 




0.09 

0.80!°;^^ 

0.066!°:°^° 


0.09 

95+° " 

"■'-'-0.38 

0.099!°;°^^ 


Ne 


24+° °^ 

"■"^^-0.04 


0.18!°;°^ 


Mg 


n 97+0.08 
"■"^'-0.06 


24+° °* 

"■"^^-0.07 


Si 


66+° '* 

"■""-0.12 


61+° °^ 

"■"^-0.12 


S 


32+° -° 

"■-'^-0.18 


31+° '^ 

"■■^^-0.16 


Ar (fixed) 


0.54 


0.54 


Ca (fixed) 
Fe 

Ni (fixed) 


0.34 

042+° °" 

"•"^^-0.009 

0.62 


0.34 
0.058!°;» 
0.62 



km s ' to allow a direct comp arison with the LETG values pub- 
lished bv lZhekov et al.l J2005h . Simultaneously we verified that 
the line fluxes change only marginally when the best-fit Doppler 
shift of 460 km s"' is used (the fluxes change by less then 5% 
which is less than the statistical errors). As can be seen from 
the LETG/RGS flux ratios presented in Table|5] all line intensi- 
ties except the O VII triplet increased over the 1 .5 years between 
the XMM-Newton (May 2003) and Chandra (Sep. 2004) obser- 



vations, consistent with the trend seen in the EPIC-pn spectra. 
Although the errors are relatively large there is an indication that 
the ionization degree of O and Ne has increased from the XMM- 
Newton to the later Chandra observations. 

As shown in Table|2]the fit with two ionization equilibrium 
models, i.e. model D, cannot really be excluded based on the 
values of x^r^ but the models A and B involving non-equilibrium 
models appear to fit the data slightly better. In either case the low 
temperature component is in equilibrium, which means that the 
ionization parameter log(n(.-t) > 12, whereas it is much lower 
by about a factor of 18 for the high temperature component 
(log(nf-t) - 10.76). This could mean that the high temperature 
component has formed predomin antly in a region with a density 
lower by that factor. We note that lZhekov et alJ (|200^ from the 
September 2004 LETG observation, which took place about 1 .5 
years later than the XMM-Newton observation discussed here, 
find the low temperature component in equilibrium and the high 
temperature component out of ionization equilibrium with an 
ionization parameter log(n(, t) = 11.23, which is a factor of about 
three higher than the XMM-Newton value. Apart from the sep- 
aration in observation time the difference between Chandra and 
XMM-Newton here is that the Chandra analysi s is restricted to 
evaluation of the emission lines and that Zhekov et a I] ll200l 
have used a different emission model, which makes use of a 
continuous distribution of shock velocities with the ionization 
parameter rigidly coupled to the post shock temperature. In the 
following we discuss a possible consequence if the high temper- 
ature component would be significantly out of ionization equi- 
librium. 

The region in and close to the equatorial plane of the inner 
ring is likely to have a higher density than the regions well above 
and below, which suggests that the low-temperature component 
is dominated by emission from the equatorial plane, whereas the 
high-temperature component is preferentially created in regions 
way off the plane. This would be an alternative to the association 
of the two temperatures with the forward and reverse shocks. A 
very simplified check can be made by comparing the emission 
measures of the two components, which show a ratio of 16/1 for 
model A. If assumed to be disk-like the fractional emission vol- 
ume of the low-temperature component scales with the ratio of 
the height of the inner ring and the radius of the inner ring. We 
take the ring height to be the same as the ring width which can 
be determined from the early optical/UV image as about 1/15 
of the radius (Panaaia & Gilmozzi 1991). Taking the density ra- 
tio as above, i.e. 18/1 and the volume ratio as 1/15 the emission 
measure ratio is expected to be 18-18/15 - 22/1 in comparison 
with the measured value of 16/1. Given the uncertainties in these 
numbers, in particular the ring width, and the spatial structure, 
which is certainly present in the matter distribution, the similar- 
ity between measurement and expectation is encouraging. 

For determining the elemental abundances we followed the 
pro cedure applied to the Chandra LETG data ( Zhekov et ^ 
2005). Because of the fairly low sensitivity of the X-ray in- 
struments values of some abundances were fixed, i.e. for He, 
C, Ar, Ca, and Ni. The abundances of the remaining seven el- 
ements N, O, Ne, Mg, Si, S and Fe are free parameters for the 
fits. Our results on the relative abundances and the results which 
Zhe kov et alJ ( l2005l) derived from the Chandra LETG measure- 
ments agree very well, despite the fact that the observations were 
taken at times differing by about 1.5 years. There are, however, 
two exceptions. We find that the abundance of Si is about 2 to 
3 times higher in the earlier XMM-Newton observations. Either 
this is related to the models or the Si abundance in the region 
of the inner ring close to its core is lower than at its inner edge. 
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Table 5. RGS line fluxes for prominent emission lines. 



Line 


Energy 


Flux' 


Ratio 




[eV] 


[photons cm"^ s"'] 


LETG/RGS 



NVII 


500 


29.0+5' X 10 


-6 


1.53±0.26 


OVII 


574/569/561 


42.9+™xI0 


-6 


0.81±0.19 


OVIII 


654 


43.5+j-=xl0 


-6 


1.17±0.13 


OVII 


666 


8.9+^JxlO- 
45.4+^-^xIO 


6 




OVIII 


775 


6 


2.26±0.69 


Ne IX 


922/915/905 


-6 


1.36±0.19 


NeX 


1022 


19.6+^:1x10 


-6 


1.86±0.32 



'Line fluxes with Itr errors. For the He-like ion triplets the sum of the 
three lines is given. 



In that sense it cannot be excluded that the entire absorption 
is due to circumstellar absorption either from a single star pro- 
genitor or a pre-explosion binary system. It is curious, though, 
that this amount of absorption is so close to the LMC total col- 
umn density. But even if the circumstellar absorption would be 
just a tenth of the total, which would be compatible with mass 
loss rates known for red supergiants, the base of the wind can 
stand close to the inner ring, which would support the sugges- 
tion that the ring is a product of the interaction of the winds of 
the red and blue supergiants. Of course structured mass loss in a 
binary is still not excluded. Whether there is such a high density, 
spherically symmetric shell just outside the inner ring radius is 
likely to be decided by observations of the soft X-ray light-curve 
over the next couple of years, when the shock wave has passed 
the inner ring. 



Besides that the XMM-Newton Si abundance (0.54 - 0.82) is 
also larger than what is generally being quoted as the general 
Si abundance in the LMC (0.31). The second mismatch occurs 
for Fe, which is found to be significantly under-abundant in both 
measurements, but the LETG derived Fe abundance is at least a 
factor of about two higher than the XMM-Newton result. Does 
the Fe abundance increase towards the core of the inner ring? 

One of the remarkable features of the inner ring is the over- 
abundance of N(1.63) and the under-abundance of 0(0.18) with 
a ratio of about 9. This ratio is consistent with the XMM-Newton 
measurements with the abundances of N ~0.8 and of O ~0.066 
yielding a ratio of about 12.1±6 (c.f. Table|5, although the in- 
dividual values of the N and O abundances are both lower by 
about a factor of 2 to 2.5. In general, the abundances of all ele- 
ments except that of He are lower than solar But it seems that 
the seven elements which were set free in the fits fall into three 
groups, which are N and Si in group one with an average abun- 
dance value of 0.73, Ne, Mg and S in group 2 with an average 
abundance of 0.28, and O and Fe in group 3 with the very low 
abundance of 0.054. Given the error boundaries there is no over- 
lap between these three groups. If this distinction is not a chance 
coincidence it might tell us something about the creation of the 
inner ring and/or the progenitor. Within each group, though, the 
abundance ratio conforms to solar values, i.e. the ratio of O/Fe is 
solar-like, as is that of N/Si, as well as the ratio for the elements 
in group 3. It is unlikely that these abundances have resulted 
from a depletion process, which leaves the creation process of 
the elements in the progenitor star as the most likely origin. 

At last we discuss the column density, the LMC Nh of 
(2.50^g33)xlO^' cm"^ for the photo-electric absorption just in 
the LMC. This can be compared with the total column den- 
sity throughout the LMC, which is measured in HI to be 
2.53x10"" cm -2 ( iBruns et alJl2005l) . The agreement means that 
either SN 1987 A is located at the far side of the LMC or that 
it shows intrinsic absorption by matter in front and close to it. 
Actually the X-ray measurement can be used to set an upper 
limit for cold X-ray absorbing matter at distances from the ex- 
plosion site greater than the radius of the optical ring, which is 
~7xl0'^ cm. If the base of an earlier wind of the progenitor star 
with a r"^ profile would be located just outside the inner ring the 
upper hmit of the base density is 4270 cm"''. For a stellar wind 
with a velocity of 10 km/s this density requires a fairly large stel- 
lar mass loss rate of 5.6x10""^ solar masses per year This value 
is larger by a factor of about six than the upper end of the mass 
loss range derived f rom radio observations of type II supernovae 
jWeileret alJ2005l) . 



4. X-ray light curve 

4.1. Data 

Since the discovery in 1994 the soft X-ray flux has been moni- 
tored first by ROSAT and then by Chandra and XMM-Newton. 
iHasinger et aL C1996.) have summarized the early ROSAT data. 
For constructing a self-consistent light curve the energy and flux 
cross-calibration of the various instruments used are needed to 
be known to a level possibly comparable to the statistical errors. 
In the following we make an attempt to re-calibrate all previ- 
ous data taken including ROSAT and Chandra. Using the best fit 
model from the EPIC-pn first-light spectra we re-determined the 
count rate to flux conversion factor to 9.05x10"^ PSPC counts 
/ (erg cm"^ s ') and also add the fluxes from five HRI and one 
PSPC observations performed up to nearly day 4000 after ex- 
plosion, which have not been published previously but are in the 
ROSAT archive. 

Chandra performed a regular monitoring of SN 1987 A since 
the launch of the observatory with published flu xes before 
the l ast XMM-Newton observation up to day 5800 jPark et alJ 
2004). Extrapolating the exponential flux increase detected by 
Chandra to the epoch of the last XMM-Newton observation we 
compute a significantly lower flux than observed with XMM- 
Newton. To investigate if this is caused by a faster brightening 
of the source or by an instrumental effect we analyzed two more 
Chandra observations which are available in the public archive 
and which were performed after the last XMM-Newton obser- 
vation. We determined count rates in the 0.5-2.0 keV band and 
apply the same flux conversion factor as iPark et alJ ll200l did 
for their last Chandra observation. The two additional flux val- 
ues follow nicely the extrapolation of the Chandra light curve 
demonstrating that the different fluxes given by the two instru- 
ments are caused by a systematic effect. Park et al. (2004) state 
in their paper that they did not correct the Chandra ACIS fluxes 
for CCD pile-up effects which they estimate to be <10%. Pile- 
up of photons causes flux losses and is therefore the most likely 
origin for the lower Chandra fluxes. In contrast the count rate 
for the EPIC-pn is far from any pile-up limit, especially during 
the SW mode observation. In addition, pile-up effects should be 
reduced for the earlier observations when the fluxes were lower, 
consistent with the better agreement between XMM-Newton and 
Chandra fluxes derived from those observations. 

We therefore suggest that the lower fluxes derived from 
Chandra ACIS data are due to pile-up and we have corrected the 
data for the flux losses following the Chandra users handbook for 
the dependence of the pile-up fraction on the photon rate with a 
normalization using the last XMM-Newton observation as refer- 
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Fig. 4. X-ray light curve of SN 1987 A 
in the 0.5-2.0 keV band. Chandra ACIS- 
S fluxes were corrected for pile-up losses 
as described in Sect. 4. The insert shows 
the light curve in logarithmic intensity 
scale together with the model described in 
Sect. 4. 



ence. We find a pile-up correction for day 6158 of about 22%, 
which can clearly not be neglected. Three more Chandra fluxes 
were published by Par ket al.l ( l2005l) . During those observations 
the CCD readout time was reduced by a factor of about two to re- 
duce pile-up eff'ects. However, for the last Chandra observations 
the intensity increased by more than a factor of two compared 
to day 6158 and therefore results in similar pile-up losses. We 
corrected also these last three Chandra fluxes as before, taking 
also into account the reduced count rate per readout frame. 

The final, corrected X-ray light curve in the 0.5-2.0 keV 
band is presented in Fig.l^with the fluxes listed in Table|6] The 
flux changes linearly with time until about day 3500 and contin- 
ues to increase exponentially until day 6700. The insert shows 
the light curve on a lin-log scale together with the prediction 
constructed from a model of the matter distribution inside the 
volume run over by the blast wave. 

So far we have concentrated on the soft X-ray supernova 
light curve but Chandra and XMM-Newton also allow to derive a 
hard X-ray light curve for the band 3-10 keV. The three XMM- 
Newton data points on days 4712, 5156 and 5918 agree very 
well with the Chandra data points given by Parket al. (2005) 
and confirm the fairly slow increase of the hard flux compared 
to the soft emission increase up to day 5918. 

4.2. Model 

The emission from SN 1987A is obviously thermal, so that 
the observed flux is proportional to the cooling function and 
the emission measure and inversely proportional to the distance 
squared. For the distance we take 50 kpc. The value of the cool- 
ing function is taken from the observational best fits to the spec- 
tra; the observed temporal changes in temperature have a very 
small effect on the cooling function in the energy band of 0.5-2 
keV, so that we take the cooling function in the band to be con- 
stant in time for the observations at hand. Therefore the shape of 
the SN 1987 A soft X-ray light curve is determined by the emis- 
sion measure, i.e. the matter density squared integrated over the 



Table 6. X-ray fluxes of SN 1987 A. 



Age 


Instrument 


0.5-2.0 keV Flux 


Pile-up 


[days] 




[10"" ergcm-2 s"'] 


correction 


1448 


ROSAT HRI 


0.091±0.113 




1645 


ROSAT PSPC 


0.204±0.057 




1872 


ROSAT PSPC 


0.250±0.057 




2258 


ROSAT PSPC 


0.363±0.068 




2408 


ROSAT PSPC 


0.431±0.091 




2715 


ROSAT HRI 


0.442±0.147 




3013 


ROSAT HRI 


0.544±0.079 




3350 


ROSAT HRI 


0.69±0.10 




3570 


ROSAT HRI 


0.74±0.10 




3640 


ROSAT PSPC 


1.01+0.09 




3650 


ROSAT HRI 


0.93±0.11 




3760 


ROSAT HRI 


0.70±0.14 




3950 


ROSAT HRI 


1.11±0.12 




4609 


Chandra ACIS 


1.68+0.06 


1.035 


4711 


Chandra ACIS 


1.81+0.07 


1.038 


4712 


XMM EPIC-pn 


1.76+0.07 




5038 


Chandra ACIS 


2.74+0.03 


1.057 


5156 


XMM EPIC-pn 


2.77+0.08 




5176 


Chandra ACIS 


3.12+0.07 


1.065 


5407 


Chandra ACIS 


4.15+0.05 


1.087 


5561 


Chandra ACIS 


4.91+0.07 


1.102 


5791 


Chandra ACIS 


6.79+0.07 


1.142 


5918 


XMM EPIC-pn 


8.10+0.09 




5980 


Chandra ACIS 


8.46+0.09 


1.177 


6158 


Chandra ACIS 


10.32+0.11 


1.216 


6359 


Chandra ACIS 


13.05+0.78 


1.136 


6533 


Chandra ACIS 


19.64+0.79 


1.205 


6716 


Chandra ACIS 


24.35+1.22 


1.255 



volume that has been encompassed by the forward shock or blast 
wave. 

The key element is to find a reasonable description for the 
spatial distribution of the matter density. We assume the in- 
ner ring to be a circular torus with a radius of r^. measured to 
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be 6.2x10'^ cm (0.83") jPanagia & Gilmozzi|[T99lt . a circular 
cross section and the origin coincident with the explosion site. 
The density distribution is independent of the azimuth along the 
torus perimeter - it is circular symmetric. Within the torus the 
density ne is described by a Gaussian, i.e. nu cc exp[-(Ar/hi)^], 
with Ar being the shortest distance from the torus center line, 
not the torus origin, though. At some torus radius r, a transition 
from the Gaussian form to a simple exponential dependence oc- 
curs, i.e. nH exp[-(Ar/ho)] for Ar > r , (see Fig.|5}. The expo- 
nential dependence is further modified by a factor of 1/r where 
r is the distance from the origin. The two forms are normalized 
to each other such that a smooth transition in density occurs at 
Ar = r,. In addition to h,, h,, and r,, which are to be determined, 
the location of the blast wave front r^ at some specified epoch 
is a free parameter. For this paper we assume the blast wave ve- 
locity to be constant in time up to day 7000. The shocked matter 
or the X-ray emitting volume elements, are assumed to expand 
radially and adiabatic with a speed of three quarters of the blast 
wave velocity. 

An acceptable fit to the data (c.f. Fig. ^ has been obtained 
for r^ = e.SxlO''^ cm (0.84") at day 7000, which means that 
the forward shock has just passed the center of the torus. This 
also means that the mean blast wave shock velocity over the 
past 7000 days is 10400 km s"' for a distance of 50 kpc. The 
transition radius r, - 3.6x10'^ cm (0.48"), so that the region 
inside the torus in which the Gaussian density distribution pre- 
vails is roughly half of the torus radius r^. The scale height is 
h; - 2.1x10'^ cm (0.28"). The scale height of the exponential 
"atmosphere" of the ring is h,, - 2.6x10'^ cm (0.35"). From 
adjusting the measured and the best-fit fluxes we get a density 
of nH,c = 1.15x10"* cm"^ for the innermost core of the ring. 
This number can be compared with the ring electron density 
of ne = (2 - 4)xl0^_cin_^ derived from the early optical/UV 
measurements bv lFransson et al. ( 1989). The ring width, which 
shows up in these optical/UV measurem ents, is 6.7x10'^ cm 
( 0.09") cPanagia & Gilmozzi 1991) and fransson et alJ jl989h 
estimate for the included mass 0.03 - 0.05 Mq. Our model gives 
for that limited region a mass of 0.065 solar masses, which is 
slightly higher. But this number is very sensitive on the exact 
width of the ring because of the Gaussian density dependence. 
The total mass overrun by the forward shock wave is, however, 
much greater and amounts to 0.45 Mq, with 90% of the mass 
contained in the much wider ring and 10% in the exponen- 
tial "atmosphere". The total mass around SN 1987 A out to day 
7000 already exceeds the previously quoted mass of the inner 
ring by a factor of 10, which may have some relevance to the 
models explaining the origin of the ring, as far as the energetics 
are concerned. If the upstream density distribution beyond day 
7000 follows the functional forms given above, a total mass of 
1.3 Mq will be involved. The soft X-ray light curve is predicted 
to flatten in mid 2006 and will reach its maximum luminosity of 
2x10^*' erg/s in early 2009. 

Unacceptable fits to the data, either in the early part or in the 
late part of the light curve, become apparent for changes of r^,, r, 
and h; exceeding 5 to 7%. This error range is so small basically 
because rc is fixed and the distribution is steep. This is different 
for ho, the scale height of the early data; h^ can be changed up 
to factor of 2 until deviations become significant. The late time 
light curve is not affected. Within this range of errors the central 
density changes by up to 30%. 

The numbers given above have been derived for a constant 
blast wave velocity v;,, which is certainly not the case according 
to the latest Chandra results (Parket al. 2004). This will affect 
our estimates of nH,c and h,-, as well as ho, because they scale 



r/ \ 
/ Ar 



Fig. 5. Sketch used to model the matter distribution around 
SN 1987 A. The matter is arranged in a circular torus of radius 
r^; lines of constant density form circles which build up the 
torus. For torus circles up to radius r, the distribution follows 
a Gaussian which has its peak at r^ and falls off with increasing 
Ar. Outside of r, the distribution is roughly exponential which 
peaks at r, and decreases with increasing Ar. For each location 
in space (r, 0, <p) the density is defined according to this scheme 
(no dependence on 0). The big circle represents the radius of the 
forward shock wave at day 7000, which has just passed the torus 
center. 



linearly proportional to v;, and inversely proportional to v^, re- 
spectively. The slow down of Vh reaches values of up to a factor 
of 2 to 3, in the late stage of the expansion, and accordingly h,- is 
expected be lower and ne.c might reach (2 - 4)xl0'*cm"^ on the 
average. In contrast, the estimates of the mass in the ring is only 
slightly affected and will stay close to the large values given. In 
a separate paper we will describe the model in more detail and 
will present results for a model with v^, slowing down in time. 

5. Limits on a compact object 

5.1. Observations 

Up to now the search for a compact remnant has not revealed any 
significant evidence for the existence of a neutron star or black 
hole, which may be due to various reasons. If a pulsar with a 
wind nebula were created in the explosion one would expect 
to see some power-law like X-ray spectrum photo-electrically 
absorbed by the ambient debris of the progenitor's envelope. 
Depending on the column density the pulsar and its nebula 
should appear at some time. To overcome the absorption one 
should look for the emission in the hard X-ray regime first, and 
XMM-Newton can supply some data and clues. 

To estimate the putative contribution from a compact object 
to the total flux a power-law component with fixed photon index 
(2.1) was added to model A in the combined fit to the EPIC- 
pn spectra. The intensity was not allowed to vary with time, i.e. 
it was constrained mainly by the first-light spectra in January 
2000. The improved only to 803.5 which formally does not 
justify the additional free parameters. For 90% confidence an 
upper limit for the flux in the power-law component can be de- 
rived. For the 0.2-10.0 keV (0.5-2.0 keV) band this results in 
an observed, i.e. not corrected for absorption, flux of 1.1x10"'^ 
erg cm"^ s"' (4.4x10"'^* erg cm"^ s"') which corresponds to an 
absorption corrected luminosity of 5.7x10-''* erg s"' (2.1x10"^'* 
erg s"') for a distance of 50 kpc. This is slightly less than the ab- 
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sorption uncorrected 2.3x10^^ erg s"' uppe r limit in the Chandra 
images (0.5-2.0 keV, ' Burrows et alJ20()Oh . Even if a pulsar was 
created in the explosion and even if the stellar debris is fairly 
transparent in the meantime there are other reasons why the pul- 
sar has not emerged in the observations. 

5.2. Discussion 

The detection of neutrinos immediately after the appearance of 
SN 1987 A is taken as proof that in this supernova a (proto-) neu- 
tron star was created. Assuming magnetic flux and angular mo- 
mentum conservation during the collapse a dipolar surface field 
of the order of lO'^ G and a rotational period not much smaller 
than 0.01 s are expected. Assuming the "standard" parameter for 
the moment of inertia / - 10^^ g cm^, we calculate the spin- 
down luminosity L^d - 4k^IP/P^ ^ 4 x 10^^ P erg s"'. Since 
the X-ray lumi nosity - 10'^ L,^ (Becker & Triimpeil l 19971 
iGil et al.'^OOl"). we find a lower limit of P > 10 in order for 
the expected X-ray luminosity of the pulsar to significantly ex- 
ceed the observed upper limit. If the rotational energy loss of 
the pulsar is due to magneto-dipole radiation the dipolar surface 
magnetic field Brf = 2x10'^ (PP-is)'^^ G, where again the "stan- 
dard" parameters have been applied, and P_i5 = P/IO"'^ s s"^. 

In order to get the X-ray luminosity of the pulsar larger 
than 5.7 x lO^'* erg s"', i.e. to enforce a spin-down faster than 
P-15 = 1, the surface dipolar field sti-ength B,/ > 2x 10'^ Vp 
G, i.e. larger than 2 x 10" G. If the field of the pulsar is lower 
it would still escape current detection capabilities. Furthermore, 
in case of SN 1987 A, matter fall -back might have been so pow- 
erful that approximately 0.1 Mq have been accreted (Chevalie3 
[1989), which would submerge the initially existing pulsar mag- 
netic fields down to regions close to the crust-core interface of 
the neutron star. Depending on the equation of state and on the 
conductive properties of the crust it takes at least 1000 years if 
not much longer, until the field is re-diffused to the sur face with 
a strength larger than 10" G rsee lGepnert eTailflggl and ref- 
erences therein). 

The only other possibility we see for the neutron star in 
SN 1987 A to have a surface magnetic field exceeding 2 x 10" 
G is the transformation of thermal into magnetic energy via the 
strong temperature gradient present in the liquid crust of the 
young star causing a thermo-magnetic instability lUrpin et'al] 
[1986; Wiebicke & Geppert 1996). The typical period to reach 
poloidal surface fields of the order of 10'^ G however, exceeds 
certainly ~ 100 yrs. In conclusion, a pulsar powered signature 
in the spectrum of SN 1987 A is not very likely to appear very 
soon. 



6. Conclusions 

We present the results from a spectral analysis of XMM-Newton 
data from SN 1987 A obtained in January 2000, April 2001 and 
May 2003. The spectra are consistent with thermal emission 
from the supernova remnant without a significant hard power- 
law component (with an absorption corrected 0. 2-10. OkeV lumi- 
nosity limit of 5.7x10-''* erg s"') originating from a compact ob- 
ject. Modeling the spectra with two components involving emis- 
sion from plasma in collisional ionization equilibrium and/or in 
non-equilibrium yields clear evidence for a temperature rise af- 
ter the second epoch observation for both components, reaching 
~0.3 keV and ~3 keV in the soft and hard component, respec- 
tively. 



We have collected all available X-ray fluxes in the 0.5- 
2.0 keV band from ROSAT, XMM-Newton and Chandra for the 
light curve of SN 1987 A up to day 6716 after the explosion. 
This includes previously unpublished ROSAT data available in 
the archive and published Chandra fluxes which we corrected for 
photon pile-up effects present in the CCD spectra. The final light 
curve shows a linear increase up to about day 4000 and an expo- 
nential rise afterwards. We model the light curve with emission 
from a circular torus around the explosion site which represents 
the inner ring around SN 1987 A. From the best fit we find that 
the forward shock just passed the center of the torus and that the 
mean shock velocity over 7000 days is 10400 km s"'. With a 
core density of the torus of 1.15x10'* cm"-' we find a total mass 
overrun by the forward shock of 0.45 M©. The model, assum- 
ing a constant blast wave velocity in its simplest form, predicts 
a flattening in the X-ray light curve for mid 2006. A new XMM- 
Newton observation is currently scheduled for the end of 2006 
to further monitor the spectral evolution of SN 1987 A. 
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